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We study ultrahigh energy astrophysical neutrinos and the contribution of tau 
neutrinos from neutrino oscillations, relative to the contribution of the other fla- 
vors. We show the effect of tau neutrino regeneration and tau energy loss as they 
propagate through the Earth. We consider a variety of neutrino fluxes, such as 
cosmogenic neutrinos and neutrinos that originate in Active Galactic Nuclei. We 
discuss signals of tau neutrinos in detectors such as IceCube, RICE and ANITA. 



1. Introduction 

Ultrahigh energy astrophysical neutrinos provide unique probes of parti- 
cle physics at energies currently not accessible to colldier experiments. In 
addition, neutrinos point back to their sources and they travel very large 
distances without interactions, providing valuable information about ex- 
treme enviroments. 

The experimental data on <-> v r neutrino oscillations 1 implies that 
astrophysical sources of muon neutrinos become sources of and v T in 
equal proportions after oscillations over astronomical distances 2 . Even 
though and v T have identical interaction cross sections at high energies, 
signals from v T — > r conversions have the potential to contribute differently 
from signals. The t lepton can decay far from the detector, regenerating 
v T 3 . In case of u^, muons produced via charged-current decay but elec- 
tromagnetic energy loss coupled with the long muon lifetime make the 
regeneration from muon decays irrelevant for high energies. Another signal 
of v T — > r is the tau decay itself 4 . 
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We have studied in detail the propagation of all flavors of neutrinos 
with very high energy (E > 10 6 GeV) as they traverse the Earth. Because 
of the high energies attenuation shadows most of the upward-going solid 
angle at high energies, so we have limited our consideration to nadir angles 
larger than 80°. We have focused on the contribution from tau neutrinos, 
produced in oscillations of extragalactic muon neutrinos as they travel large 
astrophysical distances. 

Neutrinos from astrophysical sources are usually produced via pion de- 
cays, which determine the flavor ratio v e : : v T to be 1 : 2 : 0. After 
propagation over very long distances, neutrino oscillations change this ratio 
to 1 : 1 : 1 because of the maximal <-> v T mixing. For the GZK flux, v e 
and incident fluxes are different because of the additional contributions 
from v e from neutron decay and v e from fi + decays 5 . Because of this, 
the flavor ratio at Earth is affected by the full three flavor mixing and is 
different from 1:1:1. Given fluxes at the source F® , Fjj and Fjj , the 
fluxes at Earth become: 

F Ve = Fl \ sin 2 2M2F° - F^ F°J (1) 
= F VT = \{Fl + F°J + 1 sin 2 2(9 12 (2F» - F^ F°J (2) 

where 9i2 is the mixing angle relevant for solar neutrino oscillations. We 
have assumed that 623, the mixing angle relevant for atmospheric neutrino 
oscillations, is maximal and #13 is very small, as shown by reactor experi- 
ments, as well as atmospheric and solar data. 

For GZK neutrinos 5 produced by cosmic ray interactions with the mi- 
crowave background, the flavor ratio at Earth deviates from 1:1:1 because 
the initial fluxes are somewhat different, and they start out in a ratio not 
equal to 1 : 2 : 0. In this case the v e <-» v^,v T oscillations relevant to 
solar neutrino oscillations start playing a role, in addition to the maximal 

<-> v T oscillations relevant to atmospheric neutrinos. Z burst neutrinos 6 
from models with ultrahigh energy neutrinos scattering with relic neutrinos 
to produce Z bosons are also considered below, where neutrinmixing yields 
flux ratios of 1 : 1 : 1. 

The initial fluxes for GZK and Z burst neutrinos are shown in Fig. 1. In 
our propagation of neutrinos and charged leptons through the Earth 7 , we 
have focused on kilometer size neutrino detectors, such as ICECUBE 8 and 
the Radio Ice Cerenkov Experiment (RICE) 9 and on a detector with much 
larger effective area which uses Antarctic ice as a converter, the Antarctic 
Impulsive Transient Antenna (ANITA) 10 . 
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Figure 1. Initial Neutrino Fluxes 



Signals of neutrino interactions in the rock below the ice or in ice depend 
on the energy and flavor of the neutrino. Muon neutrino charged current 
(CC) conversions to muons are noted by the Cherenkov signal of upward 
going muons in a detector such as IceCube 8 . High energy electromagnetic 
showers from v e — ► e CC interactions produce Cherenkov radiation which 
is coherent for radio wavelengths. The Radio Ice Cherenkov Experiment 
(RICE) has put limits on incident isotropic electron neutrino fluxes which 
produce downward-going electromagnetic showers 9 . The Antarctic Impul- 
sive Transient Antenna (ANITA) also uses the ice as a neutrino converter 
10 . These balloon missions will monitor the ice sheet for refracted radio 
frequency signals with an effective telescope area of 1M km 2 . All flavors 
of neutrinos produce hadronic showers. In addition, tau decays contribute 
to both electromagnetic and hadronic showers that could be detected by 
IceCube, RICE or ANITA. 

2. Neutrino Propagation 

Propagation of neutrinos and charged leptons is governed by the following 
transport equations: 

8F Vt (E, X) r °° J- NC 
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-g-a + Z^ (4) 

Here F Ut {E,X) = dN v JdE and F T {E,X) = dN T /dE are the differential 
energy spectra of tau neutrinos and taus respectively, for lepton energy E, 
at a column depth X in the medium defined by 

X = [ p{L')dL'. (5) 
J o 

For tau neutrinos, we take into account the attenuation by charged cur- 
rent interactions, the shift in energy due to neutral current interactions 
and the regeneration from tau decay. For tau leptons we consider their 
production in charged current v T interactions, their decay, as well as elec- 
tromagnetic energy loss. 

The effective decay length of produced taus does not go above 10 7 cm, 
even for E T — 10 12 GeV. This is because electromagnetic energy loss over 
that distance reduces the tau energy to about 10 8 GeV, at which point the 
tau is more likely to decay than interact clcctromagnetically 11 . 

We have found that the v T flux above 10 8 GeV resembles the flux. 
The lore that the Earth is transparent to tau neutrinos is not applicable in 
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Electromagnetic showers 
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Figure 3. Electromagnetic showers for GZK neutrinos. 



the high energy regime. Tau neutrino pilcups at small angles with respect 
to the horizon are significantly damped due to tau electromagnetic energy 
loss above E T ~ 10 8 GeV if the column depth is at least as large as the 
neutrino interaction length. 

At lower energies, E < 10 8 GeV, regeneration of v T becomes important 
for trajectories where the other flavors of neutrinos are strongly attenuated, 
but the v T regeneration is very effective. The regeneration effect depends 
strongly on the shape of the initial flux and it is larger for flatter fluxes. The 
enhancement due to regeneration also depends on the amount of material 
traversed by neutrinos and leptons, i.e. on nadir angle. For GZK neutrinos, 
we have found that the enhancement peaks between 10 6 and a fewx 10 7 GeV 
depending on trajectory. 

Fig. 2 shows the ratio of the tau neutrino flux after propagation to 
incident tau neutrino flux, for 89°, 85°and 80°. This ratio illustrates a 
combination of the regeneration of v T due to tau decay and the attenuation 
of all neutrino fluxes. For 89°, where both the total distance and the 
density are smaller, the attenuation is less dramatic, and the flux can be 
significant even at high energy. The regeneration in this case can add about 
25% corrections at energies between 10 7 and 10 8 GeV. For 85° the relative 
enhancement is around 80% and peaked at slightly lower energies, while at 
80° it is almost a factor of 3 at low energy. At 80°, however, the flux is very 
strongly attenuated for energies above a few x 10 7 GeV. It is already clear 
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Electromagnetic showers 
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Figure 4. Ratio of electromagnetic shower rates in the presence and absence of i/ M — * u T 
oscillations for GZK and l/E neutrino spectra for nadir angle 85° for a km size detector. 

from here that the total rates will be dominated by the nearly horizontal 
trajectories that go through a small amount of matter. However, rates can 
get significant enhancements at low energies where the regeneration from 
tau decays adds an important contribution even for longer trajectories. 

3. Showers 

We have translated the neutrino fluxes and tau lepton fluxes into rates for 
electromagnetic and hadronic showers at selected angles to see the effect 
of attenuation, regeneration, and the different energy dependences of the 
incident fluxes. We have focused on comparing the v T contribution to the 
v e and contributions to determine in what range, if any, z/ T 's enhance 
shower rates. Electromagnetic shower distributions for a nadir angle of 85° 
are shown in Fig. 3, while Fig. 5 shows hadronic showers. 

Fig. 4 shows the ratio of the electromagnetic shower rates at nadir angle 
85° in the presence and absence of oscillations for the GZK and Z burst 
neutrino fluxes (which have a characteristic l/E energy dependence). In 
absence of oscillations, the only contribution to electromagnetic showers 
comes from v e interactions. In the presence of — > v T oscillations, elec- 
tromagnetic decays of taus from tau neutrinos add significant contributions 
to these rates at energies below 10 s GeV. In the same time, for the GZK 
flux, v e — > v^^ oscillations reduce the number of ^ e 's at low energy, such 
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that below a few x 10 6 GeV there are fewer electromagnetic showers than 
in the absence of oscillations. 

The v T flux enhancements depend on the shape of the initial flux. The 
electromagnetic showers are more sensitive to this shape than hadronic 
ones. The relative enhancement in hadronic showers is also smaller than for 
the electromagnetic showers. This is because for the electromagnetic signal 
the only contribution in the absence of taus is from electron neutrinos, while 
for hadrons the tau contribution is compared to a much larger signal, from 
the interactions of all flavors of neutrinos. We have considered contribution 
from secondary neutrinos, which we find to be relatively small for all fluxes. 

For kilometer-sized detectors, at for example a nadir angle of 85°, the 
maximal enhancement due to v T contribution to electromagnetic shower 
rates for the GZK flux is about 50% at 3 x 10 7 GeV, while for a l/E flux, 
it is even larger, about 70%, at slightly lower energy. These energy ranges 
are relevant for IceCube. For energies relevant to RICE, tau neutrinos do 
not offer any appreciable gain in electromagnetic shower signals compared 
to v e — ■> e CC interactions, and they contribute at essentially the same level 
as to hadronic shower rates through NC interactions. 

In Fig. 6 we show the ratio of electromagnetic and hadronic shower 
rates in the presence and absence of — > v T oscillations. The maximum 
enhancement due to the presence of v T is about 40% at this angle, as ex- 
pected since this trajectory has low column density. However, as previously 
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Figure 6. Ratio of electromagnetic and hadronic shower rates in the presence and ab- 
sence of — > u T oscillations for GZK and 1/E fluxes of neutrinos for detection over 222 
km of ice. 

discussed, the enhancement occurs in a larger energy range and it peaks 
at higher energy than in the case of a small size detectors. In 1 km only 
taus with energies below a few x 10 7 GeV have a significant probability to 
decay, while much higher energy taus can decay over the total distance of 
more than 200 km. 

At very high energies tau neutrinos do not contribute large signals to 
kilometer-sized detectors because high energy tau decay lengths are very 
large, so the probability of a tau decaying in the detector is low. For 
detectors like ANITA which can sample long trajectories through the ice, 
one would expect a larger tau neutrino contribution to the signal from tau 
decay. Despite the long trajectory (222 km with a maximum depth of 1 
km for a neutrino incident at 89° nadir angle) the tau contributions to the 
electromagnetic shower rate is quite small for fluxes expected to contribute 
in the ANITA signal. For hadronic showers, the suppression of r decay 
to hadrons relative to v e NC interaction contributions is about the same 
as for electromagnetic showers compared to v e — > e. The v T contribution 
to the hadronic shower rate from interactions is about the same as the 
v e contribution. Clearly, large detectors with energy threshold below 10 8 
GeV and with very good angular resolution are needed for distinguishing 
between diffrerent neutrino flavors. 

Recently, it has been noted that rock salt formations have similar prop- 
erties to the Antarctic ice and can therefore be used as large scale neu- 
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trino detectors 14 . Salt has a higher density (p sa it — 2.2 g/cm 3 ) than ice 
(Pice = 0.9 g/cm 3 ), so it is possible to achieve an effective detection volume 
of several hundred km 3 water equivalent in salt. This is somewhat better 
than RICE but with a much smaller actual detector size. The threshold 
for detecting the radio signal from showers in salt is of the order of ~ 10 7 
GeV. Another recenly proposed experiment is is LOFAR 15 , a digital tele- 
scope array designed to detect radio Cherenkov emission in air showers. 
LOFAR has sensitivity in an energy range of ~ 10 5 — 10 11 GeV, so it can 
detect showers at much lower energies than other radio Cherenkov experi- 
ments. LOFAR will likely be configured to detect horizontal showers from 
skimming neutrinos as well. With its low energy threshold, LOFAR has 
an excellent opportunity to observe the shower enhancement at lower en- 
ergies due to v T regeneration and tau pilcup, which is not easily accessible 
in ANITA. 
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